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SECTION 1 
INTRODUCTION 


The H. E. Craaer CoaqMUty« Inc. began work under NASA Contract 
No. NAS8-31841 in March 1976 and coapleted work, including all aodlfi- 
cations to the contract, in Deceaber 1980. The work under the c<mtract 
has principally Involved the developwent of the NASA/MSFC REED computer 
code for predictii^ concentrations, dosage and deposition downwind froa 
rocket vehicle launches and the developaent and application of dispersion 
■odeling techniques for use in NASA environaental studies concerned 
with launch support and associated fuel storage and handling probleas. 

In this final report, we describe the calculation procedures and results 
of nine studies perforaed during the course of the work which were not 
included in the two formal technical reports previously submitted under 
the contract (Stephens and Duabauld, 1979; Dumbauld, Rafferty ami 
Saterlie, 1977) . 
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SEcn(»( 


SIMPLIFIED DISPERSION MODEL FOR THE 
MOMENTUM PHASE OF PLUMB EXPANSI(»I 


A Monentum Phase Dispersion Model coopute - prograa was developed 
for use on Che REEDA system for calculating concoitrations in the jet 
plines from various ducts at MSFC facilities. The asCheaaCical equations 
used in the computer program as well as the solutions for example problems 
are givoi below for the following source configurations: 

e A horizontal exhaust duct near ground level located 
at Che Marshall Hot Gas Facility (MHGF) 

• A horizontal (uchaust duct at a height of 30.5 meters 
above the ground located at the T/Acoustic Facility 

• A vertical exhaust duct at a height of 9.15 meters 
above the ground located at LRLF 

Effluent C(XDposition and other source parameters for Che three exhaust 
ducts are given in Tables 2-1 through 2-4. The concentrations at the 
duct exit of potentially hazardous gases contained in the exhaust effluent 
at each facility are given in Tables 2-1 through 2-3. Table 2-4 gives 
the remaining source parameters required for input to the models. Thres- 
hold hazard concentrations (THC) for the various constituents of the 
exhaust gases are given in Table 2-5. 

The plume expansion formula used in the calculations is based 
on expressions for Che expansion of plumes from jets described by Briggs 
(1971) where the source emissions are assumed to be conCf luous for a 
period of sixty seconds or longer. Observations of plumes from jets 
indicate that the plume radius r increases with distance from the 
exit plane according to the expression 
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TABLE 2-1 


(X»)CENTRATIONS IN THE EXHAUST 
DUCT (Xg) AT THE MHGP 













TABLE 2-2 


CONCEirrRATI(»lS IN THE EXHAUST DUCT 
(X^) AT THE T/ACCMJSTIC FACILITY 


Rtel 

Cbocentration (ppn) 

00 

BCN 

OOCI2 

HCl 

a 

OPR-488 

218.6 

32.4 

42.3 

15.6 

30.3 


209.6 

28.3 

53.9 

19.9 

38.6 

BK-250 

213.6 

30.2 

44.5 

11.3 

22.0 

FD-200 

260.3 

0 

0 

0 

0 

DC-3-6548 

267.5 

0 

0 

0 

0 

SLA-561 

294.5 

0 

0 

0 

0 


260.2 

77.1 

27.9 

10.3 

20.0 
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TABLE 2-3 


CONCENTRATIONS IN THE EXHAUST DUCT 
(Xjj) AT THE LRLF 


Riel 

Cbncentratlon (ppn) 

00 

HCN 

®Cl2 

HCl 

a 

CPR-488 

25.1 

3.7 

4.9 

1.8 

3.5 


24.2 

3.27 

4.73 

1.74 

3.39 

BX-250 

24.7 

3.5 

5.1 

1.3 

2.5 


81.4 

0 

0 

0 

0 

DC-3-6548 

83.6 

0 

0 

0 

0 

SLA-561 

92.1 

0 

0 

0 

0 


81.5 

5.0 

8.7 

3.2 

6.3 
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TABLE 2-4 

OTHER SOURCE PARAMETERS 


1 

1 

Parameter 

1 ■ - 

1 Facility 

1 _ . 

i 

msp 


T/Acoustic 

Duct Dlra^isioas 
(m) 

feigjit = 
Width = 

.363 

.629 

Diameter = 
0.762 

Hei^t = .915 
Width = .915 

Height Above 
Ground Level 
(m) 

1.32 


9.15 

30.5 

Exit Speed 
(m/sec) 

145.8 

1 


48.8 

220 

j 


Orientation of 
Exit Plane 


Hori 2 Dntal 


Vertical 


Hori 2 ontal 








TABLE 2-5 


THRESHOLD HAZARD CONCENTRATIONS FOR THE VARIOUS 
E30iAUST-GAS CONSTITUENTS 


Exhaust Gas Constituents 

Itareshold Hazard Ooncentratlcn 
(ppn) 

CD 

50 

BCN 

10 

®Cl2 

0.1 

HCl 

5 

Cl 

0.5 
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r 


X 


( 2 - 1 ) 



where 


u « ambient wind speed 

“ exit velocity at the source 
X ■ distance from the exit plane 

The radius r^ of the plime at a distance x from the exit plane is 


where 



r = effective radius of the exit 
o 

For rectangular shaped exits, the program assumes 



where 


( 2 - 2 ) 


(2-3) 


A = area of the duct exit 
o 


When the jet is directed horizontally. Equation (2-2) is strictly 
valid out to the distance at which the jet velocity becomes equal to the 
ambient wind speed. Assuming the wind direction is along the axis of the 
jet, the distance x' where v^ equals u is given by 
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-1 


(2-4) 


x’ 


3 V r 
o o 


V + 3u 
o 



If the jet is directed vertically* the plume rises due to the momentum 
generated by the jet and the stabilization height H is given by the 
expression 


where 


H 


h + 


6 V r 
o o 


u 


(2-5) 


h » height of the jet exit plane 

The distance downwind x' at which distance the stabilization height 
is achieved is given by 




+ 


2 


H 


( 2 - 6 ) 


Assuming that the material in the plume has a Gaussian distri- 
bution, the concentration at the plume centerline at height z above 
the ground and distance x<x* from the exhuast is given by 


Xjx<x'} 



(2-7) 
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Similarly* at distances x>x', the concentration is given by 


X|x>x '[ 


Xq ^o 
2 



r 


r 2 1 



r /H - z\ 


n /H + Z\ 


exp 

-.5 ] 

+ exp 

-.5 



U. I 


\ °z / 



• J 


• ^ 



(2-8) 


In Equation (2-7) and (2-8), 


X = gas concentration at the exit plane 


a_ = 




2.15 


(2-9) 


o,. = 50 a! 




X + X 


= 45 


" L 45 . 

Orix = x’f 


50 


0.9 


1/.9 


- X 


(2-10) 


( 2 - 11 ) 


al = standard deviation of the wind azimuth 
angle in radians 


X + X 


|x - x’{ 


) 


-x' 


( 2 - 12 ) 


(2-13) 


Gg = standard deviation of the wind elevation 
angle in radians 
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Concentrations at distances greater than x’ downwind from vertically 
oriented jets are calculated from Equation (2-8). The program does not 
contain provision for calculating concentrations from vertically oriented 
jets for distances less than x'. The critical distance which 

X<rac Is obtained by solving Equation (2-7) or Equation (2-8) using 
Iterative procedures to determine the value of x that Is equal to 

Requisite model inputs Include the height z and the meteoro- 
logical parameters u, and Og. To obtain maxlmxim estimates of 

for the MHGF, the height z was set equal to H (1.32 meters). In the 
remaining calculations, z was set equal to 2 meters. Concentrations 
dovnirind from horizontal jets tend to very inversely with the wind speed 
u because in Equation (2-9) increases with increasing wind speed. 
Concentrations donwwind from vertically oriented jets tend to vary 
directly with wind speed because plume rise decreases with increasing 
wind speed. For these reasons, u was set equal to 1 meter per second 
for calculations of downwind from the T/Acoustic Facility and 

MHGF and equal to 10 meters per second for the calculations of 
downwind from the LRLF. The turbulence parameters and Og were 

set equal to 0.08727 radians (5 degrees) in the calculations to take 
account of the turbulence Induced by the buildings associated with the 
three facilities. 

The values of the hazard distance calculated for the chemical 
constituents of emissions from the exhaust duct at the MHGF are given 
in Table 2-6. The largest calculated value for of 118 meters 

occurs for the phosgene (COCI 2 ) component of the ^uel texthane. It 
should be noted that no credit has been taken in the calculation for any 
increase in the height of the exhuast plume due to thermal buoyancy. 

The calculated distances may therefore be overestimates. 

Similar calculations made for the T/Acoustic Facility and LRLF 
showed that near ground-level concentrations of the exhaust gas constit- 
uents are always below the threshold hazard concentrations. 
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TABLE 2-6 


CALCULATED CRITICAL DISTAMOES FROM THE EXHAUST WJCT OF IBE MHGF 
FUR THE CHEMICAL CONSTITUENTS OF THE VARIOUS FUELS 


Riel 

Critical Distances (m) 

00 

HCN 

noci2 

HCl 

a 

All 

CPR-488 

3.7 

3.1 

103 

3.0 

30 

103 


3.6 

2.8 

118 

3.5 

36 

U8 

BX-250 

3.6 

3.0 

112 

3.2 

33 

112 

PD-200 

3.7 

- 

- 

- 

- 

3.7 

DC-3-6548 

3.7 

- 

- 

- 

- 

3.7 

SLA-561 

3.9 

- 

- 

- 

- 

3.9 


3.6 

3.3 

73 

2.0 

19 

73 

All Fuels 

3.9 

3.3 

118 

3.5 

36 

118 
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SKTION 3 

VAPORIZATION OF OXYGEN VENTED DURING LA^JNCH OPERATIONS 


Estimates of the time required for liquid oxygen to evaporate 
after venting from LOX tanks during launch operations were made for 
MSEC. The evaporation rate for liquid drops depends principally on the 
rate at which heat from the surrounding environment is absorbed by the 
drops. The rate of decrease in the drop diameter D^ caused by evap- 
oration, according to Fuchs (1959), is given by the expression 


d Do -4M D 

P ^ V V 

dt M D 

m p 

where 

= molecular weight of the evaporating vapor 

M = mean molecular weight of resulting vapor-air 
m 

mixture 

D^ = molecular diffusivity of evaporating vapor in air 
= density of air 

Pp = density of liquid 

Ap = vapor pressure of evaporating liquid at the droplet 
surface minus the vapot pressure in ambient air 


Pa Ap r 1/3 1/2-1 

— * — j2 -P 0.6 Sc . Re J (3-1) 

Pf 




"^•wrnuia) 
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” air partial pressure 
Sc » Schaldt nuBd>er 

M ~2slPa (3-2) 

Va 

« absolute viscosity of air 

Re » Reynolds number 
V D p 

. t . P a (3_3) 

Ma 

» terminal velocity of the falling drop 


To determine the time required for the drop to evaporate. Equation (3-1) must 

be integrated over all drop diameters D’ of the liquid oxygen ranging frMa 

P 

the initial diameter at venting to zero when the drop is fully evaporated. 
Since the Reynolds number and terminal velocity of the drop are functions 
of the drop diameter, the Integral Is not easily evaluated. To avoid a 
time consuming numerical integration, the Reynolds number. Re, was approxi- 
mated and taken to be constant during the evaporation period. When Re is 
held constant, the integral of Equation (3-1) yields the following expression 
for time required for the drop to evaporate 


2 
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where 


D “ average drop diameter used to determine the constant value 
P 

of Re 


We assumed the value of can be determined from the expression 


[- 


(D*)** + o’* 




A 

1 


(3-5) 


Table 5-1 presents values of the parameters in Equation (3-4) ve 
have used to solve for the evaporation time of druns with diameters of 1, 1.5 
and 3 mm. Since water drops with fiameters greater than about 3 om quickly 
break up into smaller drops and liquid oxygen has a viscosity about one- 
fifth that of water, we have assumed that the drops of liquid oxygen will 
also not be larger than about 3 mm. The results of the calculations are 
shown in Table 3-2, along with the average drop diameters from Equation (3-5) 
and the values of corresponding to these average diameters calculated 
using the McDonald (1960) technique. 
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TABLE 3-1 
INPUT PARAMETERS 


Paraneter 

Value 

N 

V 

32 g - sole ^ 

M 

m 

28.94 g • aole ^ 

D 

0.187 cm^ • s"^ 

V 

Pa 

1.2566 X 10"^ g • cm'^ 

P 

1.144 g * cm ^ 

P 

Ap 

~ 0.8 atmosphere 

Pf 

~ 1.0 atmosphere 

Pa 

1.7448 xl0^g*cm^*s^ 

T 

280.94 °K 
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TABLE 3-2 


EVAPORATION TIMES FOR SELECTED DROP SIZES OF 
LIQUID OXYGEN 


Initial Drop 

Diameter, D' 

P 

(mn) 

Average Drop 
Diameter, D 

P 

(am) 

Terminal 

Velocity, 

(cm • a S 

Evaporation 
Time, tg 

(s) 

1-0 

0.250 

97.9 

1.4 

1.5 

0.375 

160.2 

2.4 

3.0 

0.750 

320.0 

5.8 
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SECTION A 


ESTIMATES OF DWNWIMD CONCENTRATIONS FOR SPACE SHUTTLE ORBITER 
EMERGENCY LANDINGS AT WHITE SANDS MISSILE RANGE 


Hazard calculations were made for a hypothetical crash of the 
Space Shuttle Orbiter during an emergency landing at White Sands Missile 
Range. The model calculations were performed using a slightly modified 
version of a computer program developed for Dugvay Proving Ground (BJorklund 
and Dumbauld, 1977) and meteorological data based on a previous study we 
completed for White Sands Missile Range (Dumbauld and Bjorklund, 1977). 
Details of the calculations and their results are described below. 


Dispersion Model 

The peak centerline concentration downwind from an instantaneous 
source can be expressed as the product of two terms: 

Peak Centerline Concentration = {Peak Concentration Tern} 

{Vertical Term} (A-1) 

The Peak Concentration Term is defined by the expression 


KQ 


(2ir)^^^a a a 
X y 2 


(4-2) 


where 


K = parameter used to convert inputs into dimensionally 
consistent units 

Q = total source strength 

= standard deviation of the alongwind distribution of 
materials 


21 


PRECEDING PAGE BLANK NOT FH.MEO 



o * standard deviation of the crosswind distribution of 

y 

naterlal 

a « standard deviation of the vertical distribution of 
z 

material 


The Vertical Term is given by the expression 



where 


f. * fraction of the total source strength comprised of material 
in the i size category 

= fraction of material in the i^*' size category reflected at 
the ground surface (1 for complete reflection and 0 for no 
reflection) 

= depth of the surface mixing layer 

H = effective source height 
z * height above ground 

cH 

gravitational settling velocity ol the i size category 
X = distance from the source 

u » mean wind speed in the layer containing the cloud 
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For ccnveniance in writing Equation (4-3) » 0° (zero to the ze. o power) is 
defined to be equal to unity. Inspection of Equation (4-3) shows that the 
Vertical Term accounts for cloud depletion due to gravitational settling 
(V^^ > 0) and due to retention of raeterial at the surface (Y^, < 1). 


The standard deviation of the vertical distribution of ceterial 
is given by the expression 


a 


z 


oL X 
E rz 


8 


(4-4) 


where 

o' * standard deviation of the elevation wind angle in radians 
11 

= distance over which rectilinear vertical cloud expansion 
occurs downwind from an ideal point source 

8 = vertical diffusion coefficient 

X = vertical virtual distance 
z 



0_ = standard deviation of the vertical concentration distribution 
zR 

at a distance x_ downwind from the source 
Rz 
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The standard deviation of alongwind concentration distribution 
is given by the expression 


<’« ■[{&) * 


jCO 


r 


(4-6) 


where 


L{x} = alonwind cloud lenvth at distance x from the source 


0.6 Au X 




(4-7) 


; Au < 0 


Au = vertical wind-speed shear in the layer containing the cloud 




(zy - ( 2 ^)° 


(4-8) 


= standard deviation of the alongwind concentration distribution 
at the source 


u = mean wind speed at the reference height z 
R K 

p = wind profile power-law exponent 
= effective upper bound of the cloud 


H + 2.15 O ; z» < H 
z 2 m 


(4-9) 


m 


= "2 2 
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« effective lower bound of the cloud 


H - 2.15 a ; 2, > 2 

z I 


* I - 


(4-10) 


- <vH - 

; u > 

(z^ - (1+p) 

\ : S < 

It should be noted that the factor 0.6 appearing in Equatiof (4-7) is based on 
theoretical work by Tyldesly and Wallington (1965) and measurements from 
recent experiments at Dugway and other locations. 


(4-11) 



The standard deviation of the crosswind distribution of material 
is given by the expression 


a 

y 




where 


(4-12) 


o'{t} 

A 


Standard deviation of the azimuth wind angle in radians 
measured over the source emission time t 
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where 


1/5 


a;{T} 


a' {t } 
A o 


(y 


1 < T < 600 seconds 

o ~ 


(A-13) 


A o 


ry 


a 


standard deviation of the azimuth wind angle in radians in 
the surface mixing layer measured over the reference time 

distance over which rectilinear crosswind cloud expansion 
occurs downwind from the virtual point source 

crosswind diffusion coefficient 

crosswind virtual distance 


vR 

altiF ■ "^y 


; o < a'{x} X 

yR A ry 


a X 


ry 


■ 


*Ry i ”,R ^ -yy 


yR A‘‘^ ry 


Vi 


(4-14) 


- standard deviation of the crosswind concentration distribution 

ly 


at a distance x^^ _ downwind from the source 


A0' 


azimuth wind direction shear in radians within the layer con- 
taining the cloud 


tl • (lio) (^2 - ^l) 


(4-15) 


Az 


rate change of wind direction in degrees with height in the 
surface mixing layer where A0 is positive in the clockwise 
sense 
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Peak centerline time-average concentrations are obtained by multi- 
plying Equation (4-1) by the expression 




u T 

a 



(4-16) 


where 


T = 
a 


averaging time in seconds 


Source Model Inputs 

Source inputs used in the model calculations are shown in Table 

(4-1). The source strengths Q shown in the table were supplied by MSFC. 

The values of K are used to convert the source strength to cubic centimeters 

of gas at STP under the assumption that the pollutants will behave as ideal 

gases. These source strengths, when used in the models described above, will 

thus yield concentrations in units of parts per million parts of air (ppm) 

at STP. At the suggestion of MSFC, the source dimensions and 

a were calculated from the expression 
xo 



where v - total volume of gases comprising the source in cubic meters 

and under the assumption that one standard deviation of the source distribu- 
tion is equal to the radius of a sphere containing the released gases. Also 
at the suggestion of MSFC, the height of the release was set equal to 
twice the vertical source dimension. The net effect of these assumptions is 
to maximize the ground- level concentrations downwind from the source, because 
it is likely that any accident in which these fuels were released would 
result in a fire and buoyant rise of the combustior products and thus reduced 
surface concentrations. The source emission time T was set to 2.5 seconds. 

Also, because the pollutants are gaseous, the settling velocities were 

set to zero and the reflection coefficient set to unity (complete reflection). 
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TABLE 4-1 
SOURCE INPUTS 
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Meteorological Incuts 


The aeteorologicdl inputs for nighttime and daytime accidents 

shown in Table 4*-2 were selected from dispersion model inputs suggested 

for use at White Sands Missile Range (USMR) by Dumbauld and BJorklund (1977) 

and are based on measurements made at WSMR. The mixing depths H repre- 

sent rainlmua values measured at WSMR during daytime and nighttime hours 

(see Table 3-1, Dumbauld and Bjorklund, 1977). The values of a' o' and 

A E 

p shown in Table 4-2 are also based on measurements at WSMR and have 
been selected from Tables 3-4 and 3-5 of the report by Dumbauld and Bjorklund 
(1977). The nighttime values are for net radiation indices of -1 and -2 
and the daytime values axe for a net radiation index of 4-1. Thus , the 
selected inputs represent poor dispersion conditions which are intended 
to maximize ground- level concentrations. The lateral a and vertical 3 
dispersion coefficients for instantaneous sources were set to unity in the 
calculations. Finally, the wind direction shear term A6/AZ has beta set 
to zero to maximize concentrations downwind from Che source. 

Results of the Calculations 

The results of the calculations are shown in Figures 4-1 through 
4-8. Figire 4-1 shows daytime and nighttime peak centerline monomethyl- 
hydrazine concentrations downwind from the hypothetical Orbiter accident 
at WSMR and Figure 4-2 shows the peak centerline ten-minute average con- 
centrations. As might be expected, the peak centerline concentrations are 
much greater than the peak centerline ten-minute average concentrations 
near the source where the cloud passage time is shorter than Che ten-minute 
averaging time. As the cloud grows alongwind with increasing travel dis- 
tance from the source, the peak-centerline ten-minute average concentrations 
become more nearly equal Co the peak centerline concentrations. Figure 
4-3 through 4-8 show similar profiles for the remaining pollutants. 
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TABLE 4-2 


METEOROLOGICAL INPUTS 


Paraaeter 

Value 



Night 

Day 

H (•) 

30 

200 

■ 



(■ s S Zjj • 2 m) 

2 

4 

P 

.25 

.17 

(7* (radians) 

A 

.0506 

.0524 

0 ! (t * 600 s} (radians) 
A 0 

.1501 

.1571 

a 

1 

1 

e 

1 

1 

A6/AZ 

0 

0 

I 
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PEAK CENTERLINE CONCENTRATION (ppm) 



figure 4-1. Daytime and nighttime peak centerline monomethylhydrazine 
concentrations (ppm) down*’ind from a hypothetical Orbiter 
accident at White Sands Missile Range. 
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DOWNWIND DISTANCE (km) 


FIGURE 4-2. Daytime and nighttime peak centerline ten-minute average monomethyl 
hydrazine concentrations (ppm) downwind from a hypothetical Orfaiter 
accident at White Sands Missile Range. 







DAY 


10 ’ 


FIGURE 4- 



OOWNWINO DISTANCE (km) 


3. Daytime and nighttime peak centerline hydrazine concentrations 
(ppm) downwind from a hypothetical Orbiter accident at I%liite 
Sands Missile Range. 
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uhiGiNAL Page is. 
OF POOR QUAr.fTY 



PEAK CENTERLINE 

TEN- MINUTE TIME AVERAGE CONCENTRATION (ppm) 








PEAK CENTERLINE CONCENTRATION (ppm) 





a 


FIGURE 4- 


OOWNWIND DISTANCE (Km) 


3. Daytime and nighttime peak centerline ten-minute average nitrogen 
tetraoxide concentrations (ppm) downwind from a hypothetical 
Orbiter accident at White Sands Missile Range. 





PEAK CENTERLINE CONCENTRATION (ppm) 



FIGURE 4-7. Daytime and nighttime oeak centerline annnonia concentrations (ppm) 
downwind from a hypothetical Orbiter accident at White Sands 
Missile Range. 
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SECTION 5 


ESTIMATES OF DOWNWIND CX)NCENTRATIONS FOR SPACE SHUTTLE ORBITER 
EMERGENCY LANDINGS AT KENNEDY SPACE CENTER 


Hazard calculations, similar to those discribed in Section 4 
above, were also made for Orblter emergency landings at KSC using the 
dispersion model described in Section 4. 

Source Model Inputs 


Source Inputs used in the model calculations for an Orblter 
emergency landing at KSC are shown in Table 5-1. The source strengths 
Q shown in the table were supplied by MSFC. The remaining source para- 
meters were derived using the procedures outlined above in Section 4. 

Meteorological Inputs 

The meteorological inputs for three atmospheric stability cate- 
gories at KSC selected for use in this study are shown in Table 5-2. The 

values of the mixing depths H shown in the table were selected by reviewing 

m 

rawinsonde data presented in a series of reports by Stephens, Hickey and 
Greene (for example, see Stephens, Hickey and Greene, 1978). The values 
of H^ in the table are some of the lowest values found in reviewing the 
rawinsonde data for the selected atmospheric conditions and are intended 
to maximize ground-level concentrations at the longer distances from the 
accident. The values of the wind speed Uj^, power-law coefficient p, and 
standard deviations a' and a' shown in Table 5-2 also represent poor dis- 

A £< 

persion conditions for the selected atmospheric stability categories (see 
Table 5-1, Dumbauld, Bjorklund, Cramer and Record, 1970). The lateral ct 
and vertical 8 dispersion coefficients for Instantaneous sources were 
set equal to unity in the calculations. Finally, the wind direction shear 
term A0/Az was set to zero to maximize concentrations downwind from the 
accident. 
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TABLE 5-1 


SOURCE INPUTS 


Parameter 

Value 

Q (lbs) Monomethylhydrazlne 

2298 

Hydrazine 

278 

Nitrogen Tetraoxlde 

3791 

K (cm^.lb ^)Monomethylhydrazine 

2.2055 X 10^ 

Hydrazine 

3.1706 X 10^ 

Nitrogen Tetraoxlde 

1.1043 X 10^ 


6.23 

"yR <“> 

6.23 

a (m) 

6.23 

xo 


X = X (m) 

0 

ry rz 

I 

H (m) 

1 

13.4 

T (s) 

2.5 

‘‘i 

1 

Yi 

1 

-1 


\i ® > 

0 
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TABLE 5-2 

METEOROLOGICAL INPUTS 


Parameter 

Atmospheric Stability Category 

Stable 

Near Neutral 

Slightly Unstable 

H (ffi) 

100 

145 

300 

Uj^(ms at z^{2m} 

1.50 

4.00 

2.00 

P 

0.40 

0.20 

0.25 

Og(radians) 

0.0321 

0.0496 

0.0496 

o'{t =600s} radians 
A o 

0.0960 

0.1484 

0.1484 

a 

1 

1 

1 

6 

1 

1 

1 

AO/Az 

■— . 

0 

0 

0 
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Results of the Calculations 


The results of the calculations are presented in Figure 5-1 
through 5-6. Figure 5-1 shows peak centeillne concentrations of monomethyl- 
hydrajine (solid line), nitrogen tetraoxide (dashed line) and hydrazine 
(dotted line) donwwind from the hypothetical Orblter accident at KSC for 
stable atmospheric conditions. Figure 5-2 shows the peak centerline ten- 
minute average concentrations for the same pollutants and stable atmospheric 
conditions. As might be expected, the peak centerline concentrations 
are much greater than the peak centerline ten-minute average concentrations 
near the source where the cloud passage time Is shrrter than the ten-minute 
averaging time. As the cloud grows alongwlnd with Increasing travel dis- 
tance from the source, the peak centerline ten-minute average concentrations 
berome more nearly equal to the peak centerline concentrations. Figures 
5-3 and 5-4 show similar profiles for neutral atmospheric conditions and 
Figures 5-5 and 5-6 are for unstable atmospheric conditions. 
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FIGURE 5- 


Peak cencerliae concenCraCions (ppm) of monomechylhydrazine 
(solid line), nitrogen tecraoxide (dashed line) and hydrazine 
(dotted line) downwind from a hypothetical Orbiter accident at 
Kennedy Space Center in stable atmospheric conditions. 
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PEAK CENTERLINE 

TEN- MINUTE TIME AVERAGE CONCENTRATION (ppm) 



DOWNWIND DISTANCE (km) 

FIGURE 5-2. Peak centerline ten-minute average concentrations (ppm) of 

monomechylhydrazine (solid line) , nitrogen tetraoxide (dashed 
line) and hydrazine (dotted line) downwind from a hypothetical 
Orb iter accident at Kennedy Space Center in stable atmospheric 
conditions. 
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PEAK CENTERLINE CONCENTRATION (ppm) 



DOWNWIND DISTANCE (km) 


FIGURE 5-3. Peak centerline concentrations (ppm) of monoroe thy Ihydrazine 

(solid line), nitrogen cetraoxide (dashed line) and hydrazine 
(dotted line) downwind from a hypothetical Orbiter accident at 
Kennedy Space Center in neutral atmospheric conditions. 





DOWNWIND DISTANCE (km) 


FIGURE 3-4. Peak centerline cen-ninuce average concentrations (ppm) of mono- 
aethylhydrazine (solid line), nitrogen tetraoxide (dashed line) 
and hydrazine (dotted line) downwind from a hypothetical Orbiter 
accident at Kennedy Space Center in neutral atmospheric conditions. 
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PEAK CENTERLINE 

TEN-MINUTE TIME AVERAGE CONCENTRATION (ppm) 



DOWNWIND DISTANCE (Km) 

FIGURE 5-6. Peak centerline ten-minute average concentrations (ppm) of mono- 
nethylhydrazine (solid line), nitrogen tetraoxide (dashed line) 
and hydrazine (dotted line downwind from a hypothetical Orblter 
accident at Kem.edy Space Center in unstable atmospheric condi- 
tions. 
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SECTION 6 


HCl CONCENTRATIONS DOtfMUIND FROM ACCIDENTAL SRM IGNITION 
IN THE VAB AT KENNEDY SPACE CENTER 


A preliffllnary analysis of HCl concentrations dowwind fron the 
VAB has been sade for the foUovlog four types of postulated accidents: 

e Type I — Ignition of a single, plugged, center seg^ 

B^t in High Bay 4 with the High Bay doors closed 

e Type II — Ignition of a single, plugged, coiter 
segment in High Bay 4 with the High Bay doors 50 
percmt open 

• Type III — Ignition of a three stack (2 caiter, 1 
aft) group (unpli^^) on the MLP in High Bay 3 with 
the High Bay doors closed 

• Type IV — Ignition of a three stack (2 center, 1 aft) 
group (unplugged) on the MLP in High Bay 3 with the 
High Bay doors 50 percent open 

The source information for these four types of accidents was 
calculated for Dr. Koller by Dr. Beniss, PRC System Services Company. 

Table 6-1 shows the parameters derived by Dr. Beraiss that are required for 
calculating HCl concentrations downwind from the VAB. 

Cloud Rise Calculations 


According to the guidelines published by the ASME (Smith, 1968) 
a cloud rise exceeding 1.5 times the building height would result in a cloud 
stabilization height above the region of disturbed flow and thus greatly 
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TABLE 6-1 


S(Mnt£E PARAMEIQtS FOE THE POSTDLATED 
TYPES OF ACCinOITS 


Source 

Accident Type 

Paraaeter 

I 

II 

III 

IV 

Gas Exit ^ 

Teaperature ( C) 

167 

167 

458 

458 

Gas Exit Voliaie 

(■3) 

2.1520 

3. 1994 

7.3346 

9.3497 

Gas Exit Area 
(.2) 

612 

1955.4 

612 

1955.4 

Release Tiae 
(s) 

932 

932 

847 

847 

Exit Velocity 
(a s-h 

3-77 

1.76 

14.15 

5.65 

HCl Release 
Rate (kg s~^) 

10.172 

15.12 

61.596 

78. 144 


50 










reduce ground-level concottration levels close to the source. For this 
reason* we have used the source Inforaatlon in Table 6-1 to calculate 
cloud rise for the four types of accidents outline above. 


For the purposes of this study* cloud rise is calculated from 
the following expression (Briggs* 1971) for ataospheres with {Mtential 
temperature gradients greater than zero: 



/Sh •A 




f 


3F 

u{h}Y^s 




uu{h}s“^^^ < lOh 


TTuftJs"^^^ > lOh 


(6-1) 


where 


Ah > cloud rise 



g = acceleration due to gravity 
A * exit area of the gas cloud 
w = gas exit velocity 

» ambient air temperature 
* gas exit temperature 


( 6 - 2 ) 
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A 36 

T 3z 


(6-3) 


36 

“ vertical poteatial teaperature gradient 


entrainnent coefficient 


correction factor linitlng cloud rise as the mean wind 
speed at stack height approaches or exceeds the exit 
velocity w 


u{h} 




< W 


; u{h} > w 


(6-4) 


u{h} 


mean wind speed at the building height h 


u{2m}(|) 


(6-5) 


u{2m} ” mean wind speed at a height of 2 meters in an area not 

influ^ced by the building 


p ” wind power-law coefficient 


The source and metelrological parameters required to calculate 
cloud rise using Equation 6-1 are respectively shown in Tables 6-2 and 
6-3. The area A shown in Table 6-2 is for 4 roof ventilators commonly 
found in a grouo on the roof of the VAB. Inspection of Equation 6-1 
shows that the cloud rise increases as the value of A Increases, if all 
other parameters are held constant. Since we wished to tsake conservative 
estimates of cloud rise, we selected the area of 4 ventilators as a minimum 
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TABLE 6-2 


SOURCE PARAMETERS FOR CALCULATING 
CLOUD RISE 


Source 

Accident Type 

Parameter 

Type I 

Type II 

Type III 

Type IV 

A imh 

23.5 

23.5 

23.5 

23.5 

w (as”^) 

3.77 

1.76 

14.15 

5.65 

h (m) 

160 

160 

160 

160 

T 

s 

440.2 

440.2 

731 

731 

Y 

C.66 

0.66 

0.66 

0.66 


TABLE 6-3 

METEOROLOGICAL PARAMETERS FOR CALCULATING 
COULD RISE 


Meteorological 


Wind Speeds 
•j{2m} (ms~^ 

) 

Parameter 

1 


4 


(Very Light) 

imgii 

(Moderate) 

T^(°K) 

298.2 

298.2 

298.2 

(deg m'S 
oz 

0.02 

0.01 

0.005 


0.25 

0.20 

0.15 
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area through which the exit velocity w vas equivalent to the value shown 
In the table. The value of 0.66 for Y Is the value suggested by Briggs 
(1972) for stacks. We made cloud^rlse calculations for three nominal 2-meter 
wind speeds* as shown In Table 6-3. Values of 30/6x and p were assigned 
to these mean wind speeds on the basis of previous experience and knowledge 
of meteorological tower measurements made at KSC. 

The results of the cloud-rise calculations are given In T{d>le 6-4 
and show zero cloud rise Is calculated for Type II accidents* for Type I 
accidents with mean wind speeds of 2 and 4 meters per second and for Type 
IV accidents with a mean wind speed of 4 meters per second. However* 
none of the cloud rise values calculated for the other cases results in 
a cloud stabllizatlcn height greater than 2.5 times the building height 
(400 m) which Is the upper boundary of the region of disturbed flow as 
specified in the ASME guidelines. It should be noted that if we assumea 
the exit velocity to apply to a larger cross-sectional area or considered 
the enhancement of plume rise due to the combining of plumes from other 
groups of vents (see Briggs* 1974)* we could have calculated stabilization 
heights greater than 400 meters for some cases. However, because of 
the uncertainties associated with these plume-rise calculations and to 
preclude overestimates of the hazard distances* we have included building 
wake effects in all the model calculations. 

Building Wake Effects 

To our knowledge, the best quantitative description of building 
wake effects on concentration patterns downwind from buildings is provided 
by wind-tunnel tests conducted by EPA (Thompson and Lombardi, 1977). The 

results of the EPA study are very similar to results reported by Robins 
and Castro (1977a, 1977b). The VAB building dimensions provided by 

Dr. Roller (h = 160m, w - 156ra, d » 127m) show that the VAB building 
is approximately cubical. Therefore, we have use of the results presented 
by Thompson and Lombardi for a cubical building to evaluate the wake effects 
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TABLE 6-4 


CALCULATED CLOUD RISE (Ah) AND CLOUD STABILIZATION HEIGHT (h + Ah) 
IN METERS FOR FOUR TYPES OF POSTULATED ACCIDENTS 


Mean 2-meter 
^Wind Speed 
u{2m} (ms~^) 

Accident Type 

I 

II 

III 

IV 

Ah 

h+Ah 

Ah 

h+Ah 

Ah 

h+Ah 

Ah 

h+Ah 

1 

53.3 

213.3 

1 

160 

162.7 

322.7 

119.8 

279.8 

2 

0 

160 

B 

160 

175.0 

335.0 

57.9 

217.9 

4 

0 

160 

B 

160 

188.2 

348.2 

0 

160 
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of the VAB building with the wind direction at a 45-degree angle to a build- 
ing wall (their Case 2). This wind direction produce > the highest ground- 
level concentrations downwind froa a cubical building. 


As explained below, we used the data for Case 2 in the Thompson 
and Lombardi report to develop the model source inputs and the lateral 
and vertical dispersion coefficients to be used in our model calculations. 
The wind-tunnel measurements of vertical concentration profiles for Case 
2 (see Figure 12 on page 22 of the Thompson and Lombardi report) show that 
the plume centerline is brought to the ground by the building wake circu- 
lation at a distance of three to five building heights. To simplify our 
model calculations, we assumed a ground- level source located at the base 
of the VAB building with an initial vertical dimension (2. 15 a _) given 
by the height of the VAB building and an initial lateral dimension (-t . 30 
given by the diagonal of the horizontal cross section of the VAB 
building. Values for the lateral and vertical turbulent Intensities to 
be used in the model calculations were selected by a cut-and-try procedure 
where the objective was to obtain a normalized ground-level concentration 
profile from the model calculations, using the above source parameters in 
our quasi-continuous dispersion model described below, that matched the 
non .lized ground-level concentration profile for Case 2 at distances 
from 5 to 20 building heights. Using the source parameters described 
above, with the lateral and vertical dispersion coefficients thus deter- 
mined in the quasi-continuo'<s dispersion model, we were also able to cal- 
culate vertical concentration profiles at distances of 5 and 8 building 
I eights and a lateral concentration profile at a distance of 5 building 
heights which were in good agreement with the corresponding measurements 
presented by Thompson and Lombardi (see Figures 1? and 13 in their report) 
for Case 2. 


asi- Continuous Source Model 

The Cuasslan model for a quasi-continuous source emitting it a 
constant rate over the time t is given by the expression 

D 
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|5K-'C"F)1-[im’I 


Q = source eaission rate 


- st.-uidard deviation of the crossvind concentration distribu- 


(x + X - X (l-aj 


T jr 

A ry 


- standard deviation of the wind azimuth angle in radians 


X - downwind distance from the source 


X =• crosswi-.id virtual distance 

y 
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X, + a - a) 


(»-«) 


X 

y 


ax 

ry 



l/a 


a « crossvind dispersion coefficient 


X " distance over which rectilinear crossvind cloud expan- 

sion occurs downwind iron a virtual point source 


o <■ standard deviaticm of the crosswind concentration dis- 
~ tribution at a distance x^ fora the source 

o “ standard deviation of the vertical concentration distrib- 
ution 



standard deviation of the wind elevation angle in 
radians 


x^ * vertical virtual distance 



• stiuidard deviation of the vertical concentration dis- 
tribution at x„ 

K 

u mean wind speed 


( 6 - 9 ) 


( 6 - 10 ) 
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« mean wind speed at the referoice hei^t 
p • wind profile power-law exponent 
Zj » effective lomr bound of the clowl 



Z2 * effective upper bound of the cloud 


[H + 2. 15 a ; z- < H 
z z m 


n 


; 2-. > H 


y * crosswind distmce from the cloud centerline 
H = depth of the surface mixing layer 
H * effective height of the source above ground level 
« h + Ah 


( 6 - 11 ) 


( 6 - 12 ) 


(6-13) 


(6-14) 
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h ” actual source height above ground level 
Ah * cloud rise 
z “ height above ground level 

t “ tiae after release begins at the top of the roof 
tg ■ source eaission tiae 

o “ standard deviation of the aloogwind concentration dis- 
tribution 



Lfx^ 


0.6 A u 

u 

0 


X 


Au ^ 0 
Au < 0 


Au = vertical wind-speed shear in the layer containing the 
cloud 



o 


xo 


standard deviation of the alongwlnd cloud distribution 
at the source 


Source and Meteorological Inputs for 
the Concentration Calculations 


THr. r., ‘ concentration 

model given by Equation (6-o) are presented in Tables 6-5 and 6-6. The 


(6-14) 


(6-15) 


(6-16) 
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TABLE 6-5 

S(HJRCE INPUTS FtHt THE ONiCENTRATlON CALCULATIONS 


Source 


Accident Type 


Paraaeter 

I 

II 

III 

IV 

3 -1 

Q(ppa m s 

6.821x10® 

1.014x10^ 

4.131x10^ 

5.240x10^ 


932 

932 

847 

847 

H (o) 

0 

0 

n 

0 

V”«o <■> 

46.9 

46.9 

46.9 

46.9 


74.4 

74.4 

74.4 

74.4 


0 

0 

0 

0 


TABLE b-6 

METEOROLOGICAL INPUTS FOR THE CONCEZfTRATIOH CALCULATKNSS 


Jfeteoro logical 


Wind Speed 


Paraaeter 

Very Light 

Light 

Moderate 

o' (radians) 
A 

0.1309 

0.1309 

0. 1309 

Og (radians) 

0.0400 

0.0400 

0.0400 

a 

0.9 

0.9 

0.9 

Ujj{Zj^*2ra} (a s”^) 

1 

2 

4 

P 

0.25 

0.20 

0.15 

m 

200 

400 

800 
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source ealssloa rates In Table 6-5 vere obtained by assiaiing that 1 kilogram 
per cubic meter of HCl at a pressure of 1013.2 millibars and te^>erature 
of 25^ Celsius is equivalent to 6.7058 x 10^ parts of BCl per million parts 
of air and then multiplying by the HCl emission rates specified by KSC 
(see Table 6-1). The source emission times t^ were also specified by KSC. 

As noted earlier, we determined from an analysis of the mea- 
surements presented by Thompson and Load>ardl (1977) for a cid>ical building 
oriented at a 45-degree angle to the mean wind direction that the normal- 
ized c<mcentration profile at distances beyond 5 building heists was best 
represented by assuming a virtual volume source located at ground-level 
(8*0). This assumption was also found by Robins and Castro (1977b) to 
apply to their wind-timnel concentration measur^rats. The lateral and 
alonguind source dimensions o _ aisi o were both set equal to the diag- 

XO 

<mal length of the building (201 neters) divided by 4.3 umier the assump- 
tion that the crosswiod aid almigwind ccmcentration distributions at the 
source were Gaussian. Because the source is located at ground-level, the 
vertical source dimension was set equal to the building height h 
divided by 2.15 which follows from the assumption that the vertical con- 
cent rat ion distribution is also Gaussian. Finally because the source 
dimaision were defiled at the downwind base of the building, the reference 
distance was set to zero. 

Table 6-6 shows the values for the lateral dispersion coefficient 
0 ^ and the vertical dispersion coefficient 0^ obtained from the cut-and- 
try procedure using the dispersion model given by Equation 6-6 with the 
above source parameters, to calculate normalized concentrations equal to 
the normalized concentrations given by Thompson and Lombardi for Case 2 at 
downwind dis*‘''nces of 5 to 20 building heights. Also, the 2-meter wind 
speeds and ..le power-law coefficients in Table 6-6 are the same as those 
used in the cloud-rise calculations described earlier. Values for the 
depth of the surface mixing layer shown in Table 6-6 are estimates 

we believe to be representative of conditions at KSC for the three mean 
wind speeds. 


62 



Results of the Calculations 


Figures 6-1 throi^h 6-4 show the calculated grow^-level concen- 
trations downwind froa the VAB for the four types of postulated accideats 
dtiring very light, light and Moderate wind speeds. 

The results of the conc«tration calculations are sunBarized in 
Table 6-7. The «axl«B HCl c<»centration occurs in every case at 275 to 
300 Meters downwind from the VAB. Inspection of Table 6-7 shows that 
ground-level concentrations are greatest, as expected, for Type IV acci- 
dents. Under very light wind spcieds, HCl concentraticms above 100 ppm 
extend to a distaiu:e of about 4 kilooeters downwind froa the VAB ai^ caa- 
centratlons above S ppa extend to a distance of about 18 klloaeters down- 
wind froa the VAB. 

He also calculated tiae concentration profiles for those cases 
in which the HCl concentration exceeded 100 ppm at distances greater than 
800 Meters froa the VAB. In each of these cases, the tine profile was 
calculated at the approximate distance where the HCl concentration was 
equal to 100 ppm. The results are shown in Figure 6-5 through 6-10. 

Figure 6-5 shows, for example, the tiae concentration profile at a distance 
of 3000 meters froa the VAB for a Type III accident under very light wind 
speeds. At this distance, the HCl concentration exeeeds 5 ppa at 16.2 
minutes after the cloud starts to exit the root vents and remains above 
5 ppm until 38.2 minutes after the cloud starts to exit the roof vents. 

Figure 6-11 illustrates the effect on the calculated ground-level 
HCl concentrations of including plume rise and neglecting building wake 
effects for a Type IV accident under light wind speeds. As shown in the 
figure, the naxinua HCl concentration is reduced in magnitude from 1046 
ppm to 40 ppm and the distance to the maximum is increased from about 275 
meters to 3.3 kilometers if building wake effects are neglected. The thres- 
hold distance to an HCl .oncentration of 5 ppa is app' oximately 20 kilo- 
meters in both cases. 
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CONCENTRATION (ppm) 



DOWNWIND DISTANCE (km) 

FIGURE 6-1. Ground-level HCl concentrations downwind from the VAB for a Type 
I postulated accident. 


. page 
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CONCENTRE 






CONCENTRATION (ppm) 



DOWNWIND DISTANCE (km) 

FIGURE 6-4.. Ground-level HCl concentrations downwind from the VAB for a Type 
IV postulated accident. 
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TABLE 6-7 


MAXIMUM HCl COHC: ^RATIONS AHD THRESHOLD DISTANCES TO HCl 
CONCENTRATIONS OF 100 and 5 PARTS PER MILLION (ppit) 


Accident 

Type 

’ -1 - ■ . . , 

Hind 

Speeds 

Maximum 

Concentration 

(ppm) 

Threshold Distance 
(meters) 

100 ppm 

5 ppm 

I 

Very Light 

231 

575 

9200 


Light 

136 

400 

3900 


Moderate 

83 

* 

2900 

II 

Very Light 

344 

750 

11500 


Light 

202 

510 

5200 


Moderate 

123 

360 

3600 

III 

Very Light 

1401 

3000 

25000 


Light 

824 

1650 

16500 


Moderate 

501 

1070 

8500 

IV 

Very Light 

1777 

3950 

28000 


Light 

1046 

1950 

19000 


Moderate 

635 

1320 

9600 
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lO meters 










FIGURE 6-6. Time profile of HCl ground-level concentration at 1650 meters downwind from the VAB for 
a Type III postulated accident during light winds. 
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FIGURE fi- 
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FIfJURE 6-8. Time profile of HCl ground-level concentration at 3950 meters downwind from Che VAB for 
a Type IV postulated accident during very light winds. 
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FIGURK 6-9, Time profile of HCl ground-level concentration at 1950 metcra dovmwlnd from the VAB for 
a Type IV postulated accident during light winds. 
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CONCENTRATION (ppm) 



DOWNWIND DISTANCE (km) 

FIGURE 6-11. Ground-level HCl concentrations downwind from the VAB for a Type 
IV postulated accident during light winds. 
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SECnOH 7 


CALCUIATIOHS OF HYDRAZINE aWCENTRATIONS DOWNWIND FROM THE 
ORB ITER PROCESSING FACILITY HYPER EXHAUST STACK 


Calculatloos of hazard distances for the release of hydrazine 
from the Orbiter Processing Facility (OFF) Hyper Exhaust Stack were nade 
using the quasi-continuous source aodel described above in Section 6. The 
exhaust stacL has a height of 18.3 aeters and a diaaeter of 1.83 aeters. 

The flow rate through the stack, supplied to us by MSFC, is 434.2 a^ s 
at a teaperature of 23.89^C. According to MSFC, 1 liter of hydrazine is 
assiffled to be spilled. Since no evaporation rate was supplied, we have 
nade the assuaption that the 1 liter of hydraziitt is evaporated in 1 secoiui. 


Since the stack gases are near ambient temperature, only aomentum 
forces contribute significantly to pluae rise. Under this assumption, the 
final rise of the plume (see Dumbauld, Bjorklund and Bowers, 1973, p. A-li, 
Equation A-44) is 


h 


+ 





r 

o 



(7-1) 


where 

w = stack exit velocity = 16.53 a s 

o 


r ^ stack exit radius * 0.9144 a 

o 

Y * entrainment coefficient = 0.6 

c 


s = stability parameter 

= £ M 

T Az 


( 7 - 2 ) 
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hot filmed 



g “ 9.f ai 

T > absolute teoperature • 297.05^ 

■ la|»e rate of potential temperature 
At 

u « aean wind speed 

h > stack height • 18.3 ■ 

For worst case cc»ulitl<ns, we have set A9 /Az equal to 0.01S°K ■ ^ and u 

equal to 2 a s Thus, the pluae froa tl» OFF stadc stabilizes at a 
height of 51.74 aeters and has a radius r equal to 


r * - h) + r ” 20.99 a 

o 


(7-3) 


As noted above, the quasi-continuous source aodel described in 

Section 6 was used to calculate concentrations downwind froa the stack. 

The source Inputs based on the ^K>ve assuaptions and required for the 

calculations are given in Table 7-1. The source diaensions o a and 

2K yK 

were calculated fr<» the expression 


zR 


yR 


xo 


r 

4.3 


20.99 

4.3 


9.76 


(7-4) 


The aeteorological inputs used in the aodel calculations are shown in Table 
7-2. A nominal value for of 300 meters was used ir civ- calculations. 

The value of the wind power-law coefficient p was set tj zero, which means 
that the wind speed of 2 a s ^ is invariant with he c under the worst- 
case conditions. Because of the proximity of the OFF stark to the VAB, 
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TABLE 7-1 


SOURCE INPUTS 


Paraaeter 

Value 

Q (ca^ s“^) 

1000 

K 

1 

°.ii <■> 

9.76 

°yR <■> 

9.76 

<■> 

9.76 

•sy ■ *Rz <"> 

0 

*iy ■ Sz <•> 

50 

H (m) 

51.74 

T (s) 

1 
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TABLE 7-2 

METE(»M)LOGICAL INPDTS 


P«raaecer 

Value 

(■) 

300 


2 

P 

0 

(radians) 

0.0971 

o' {x «600s} (radians) 

A O 

0.3A91 

a 

1 

0 

1 


0 

hz 
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the value of was set equal to 20 degrees (.3491 radians) and the 
value of Og was set equal to o^{t«1 s), or 5.56 degrees (.0971 radians) 
to account for the mhanced turbulence caused by the VAB. 

The results of the calculate ns are presented in Figure 7-1. The 
highest concentration of 1.7 x 10 ^ ppn occurs at a distance of 300 
neters downwind fron the stack. This concentration is well below the 
critical level of 3 ppa for public and occupational emergency expo- 

sure limits. It should be noted that the concentration of N.H, in the 

2 4 

stack. 


3 

1000 cm N-H, per second 

2.3 ppm, 

434.2 a per second 

is also below the critical exposure limits even when the 1 liter of ^ 2 !!^ 
is released '>ver a 1-second time period. Thus, these calculations show 
that the release of 1 liter of 1 * 2^4 through the OFF stack does not 
represent an exposure hazard. 
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GROUN 



FIGURE 7-1. Preliminary estimates of concentrations at ground-level down- 

wind from the OFF Hyper ExHaust Stack. 
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SECTION 8 


PEAK AND TIME-AVERAGE HCl CONCENTRATIONS FROM 
ACCIDENTAL IGNITIONS OF SRM SEGMENTS 
IN STORAGE AREAS 

Peak and time-average HCl concentrations resulting from accidental 
Ignitions of SRM segments stored in open areas at Kennedy Space Center (KSC) 
have been calculated using a quasi-continuous source model for time-average 
concentrations. 

The calculations were performed for various source configurations 
specified by Dr. Koller at KSC according to the number of segments involved 
in the conflagration. The total solid fuel weight of the 4 segments com- 
prising a single SRM are shown in Table 8-1. Each segment is 12 feet in 
diameter. In their storage configuration, the 4 segments for SRM can 
be stored in a line with a separation distance (center-to-center) of 34 feet. 
When SRM's are stored in parallel lines with a similar spacing between 
segments and a separation distance (center-to-center) of 31 feet between the 
two parallel lines. The burn configurations we have considered in the 
calculations are shown in Table 8-2. In addition to the SRM and 2 SRM configu- 
rations, we have considered 1 forward center SRM stored alone and another 
configuration in which the 2 aft segments are stored separately from the 
other 6 segments when 2 SRM's are stored. In the latter case, either the 
2 aft segments or 6 other segments are assiuned to burn separately. 

Meteorological Conditions 

The concentration calculations were made for very light "worst- 
case" wind conditions and for normal wind conditions at KSC. The meteorologi- 
cal conditions used in the plume-rise and concentration models described 

below are shown in Table 8-3. The model parameters, except for H , shown in 

m 

the table for very light wind conditions are representative of nighttime 

stable conditions at KSC. The depth of the surface mixing layer H for 

m 

these conditions is normally lower than the value of 550 meters shown in 
the table. However, if we had set H to a T jwer value, the buoyant plume 
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TABLE 8-1 

SOLID FUEL WEIGHT FOR SRM SEGMENTS 


Segment 

Fuel Weight (lbs) 

Forward 

301, 125 

Forward Center 

273, 318 

Aft Center 

271, 977 

Aft 

261, 277 


TABLE 8-2 

BURN CONFIGURATIONS USED IN THE CALCULATIONS 


Burn Configuration 

Segments 

1 

Forward Center 

2 

Two aft segments 

4 (1 SRM) 

One each forward, forward center. 


aft center and aft segments 

6 

Two each forward, forward center 


and aft center segments 

8 (2 SRM’s) 

Two each forward, forward center 


aft center and aft segment" 
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TABLE 8-3 

METEOROLOGICAL MODEL INPUT ’ A •A>iETERS 


Model 

Wind Conditions 

Parameters 

Very Light 

r 

Moderate 

u{2m} (m s 

1 

4 

P 

0.25 

0.15 

(m) 

550 

800 

0.{x *10 min} (deg) 
A o 

5 

8 

(Jg (deg) 

1.7 

2.7 

A$/Az (°K 

0.02 

0.005 


0.66 

0.66 

T (°K) 

298.2 

298.2 

Pa 

1183.9 

1183.9 
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froa soae of the burn conf Jgurai:loa8 would have peMtrated the elevatml 
inversion at the top of the aizing layer and the pluae would not have aixed 
to the ground. Because the objective of this stwly was to deteraine the 
aaxiaua hazard distance under "worst-case" conditions, c)m value of used 
in the calculations was selected to yield a aezizaia hazard dis*-ance. The 
Bodel parameters shown in Table 8-3 for moderate wind conditions represent 
noneal daytiae conditions at KSC. Definitions of the input parameters are 
given in the discussions of the plume-rise and dispersion models below. 

Pluma-Rise Calculation^ 


Pluae rise for the quasi-continuous sources forMd by the Iwming 
segnents was calculatctd from a model developed by Briggs (1970) for contimiCHis 
source emissions into a stable atmosphere- The derivation of the pluate-rise 
model is givcm in Apper:'ix A of KASA CR-129006 (Du^auld, Bjorfclund and 
Bowers, 1973). The height z of the cloud centroid at time t after release 
is given by the expression 

r 3 1 


Ah 



s 




( 8 - 1 ) 


where 

F = 
c 


gQc 




( 8 - 2 ) 


- 2 . 

g = gravitational acceleration (9.8 a s “) 


= effective rate of heat release (cal s S 

-3 

* air density (g n ) 


Cp = speci‘"ic heat of air ac conctant pressure 
(.24 cal 
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T “ air cenperature (°K) 

s “ stability paraa«ter 

. & ^ 

T AZ 


(8-3) 


^ « vertical gr^ient of ambient potential te^[>erature 

u “ mean wiiui speed fr(» a reference height of 2 meters to 
the height z 


(z-2) 2** (1+p) 
radius nf t!^ source 
extraxment parameter (0.66) 
The final plume rise from Equation (8-i) is 



%diich occurs at the time 



(S-4) 


(8-5) 


t 



( 8 - 6 ) 


Inspection of Equation (8-1) and i8— i) shows that the height z is dependent on the 
mean wind u which is in turn dependent on z. Thus, iteration is required to 
determine the plume rise at any time t. 
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Use of t}ie plvae-rlse aodel requires that aod r^^ be specified. 

Values of Q were obtained by assuainx that the burning fuel releases 1100 
c 

BTU's per pound, or ab<Hit SO percent of the stoichioaetric ualue of heat 
available. Also, each segaent was assuaed to bum at a constant rate over 
a period of 932 seconds. Table 0>4 shows the effective value of used 
in the pluae-rise calculations for the various bum configurations. The 
values for the initial pluee radious r_ shown in Table 8*'4 were calculated 
under the assuaption that the area A containing the burning segaents for 
each configuration could, for our purposes, be represented by a circle 
with an equivalent area aod radius r^^ given by the expression 

(8-7) 

Finally, Table 8-4 shows the total aaount of BCl released over tiw: 932-second 
period which was based <sa ti» assuaption that 21 percent of the fuel weight 
in each segaents was released ^ gaseous HCl. 

Tables 8-5 and 8-6 show the results of using the plwe-rise aodel 
to calculate the centerliiK height of the pluae and the pliae radius at 
that height (r > Y ^h^tg) at various downwind distances froa the Immlng 
segaent(s) for, respectively, very light aiul aoderate wind conditions. 


Quasi-Continuous Source Model for 
Tine-Average Concentrations 


Tne 

downwind from 
emission time 

Xtx, 


Gaussian aodel for calculating tiae-average concentrations 
a quasi-continuous source with a constant rate over the source 
Tg is given by the expression 


■> 


Ta+b* 

Q 1- 

\ 1 / y 1 1 1 

1 2x 

- P n 


2tt 0 a u » 

w •* 



y * 


( 8 - 8 ) 


(Equation (8-8) continued on page 88) 
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TABLE 8-4 
SOURCE PARAMETERS 


Segment 

i 

Effective Heat i 

Release Rate, ! 

(g cal s 

— 

Initial 

Plwe Radius, r_ 

(■) 

RCl Enlssion 
Rate, Q 

(r s“^) 

1 

8.1288 X lo' 

3.7 

2.7935 X 10* 

2 

1.5541 X 10* 

4.9 

5.3408 X 10* 

1 4 (ISBM) 
1 

3.2943 X 10® 

9.5 

1.1321 X 10^ 

! ^ 

5.C345 X 10* 

12.2 

1.7302 X 10^ 

1 8 (2S»!'s) 

t 

1 

\ 

6.5886 X 10* 

14.3 

2.2643 X 10^ 
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TABLE 8-5 

PLUME USE PABAMETEBS FOR VERY LIGHT tflHD OWDITICRIS 


Downwind 

Distance 

Centerline Height (a) 


Plume Radius 

(•) 


Hiad>er of Segments 

Ihodier of Segments j 

(■) 

l_lJ 

2 

4 

6 

8 

1 

2 

4 

6 

8 

100 


217e7 

255.3 

284.3 

304.3 

121.2 

143.7 

177.9 

199.9 

215.2 

125 


242.8 

286.2 ’319.1 

1 

341.9 

135.0 

160.2 

198.3 

222.9 

239.9 

160 


272.1 

322.6 

360.3 

386.4 

150.9 

176.6 

222.4 

250.1 

269.4 

i 200 

244.8 

298.8 

356.3 

398.7 

428.1 

165.2 

197.2 

244.6 

275.4 

296.9 

1 250 


323.5 

388.1 

435.5 ^ 

468.4 

177.9 

213.5 1 265.6 

299.7 

323.4 

i 320 

277.1 

342 6 

415.2 

468.2 

505.0 

186.5 

226.1 283.5 

1 

321.2 

347.7 

, 338.2 

277.7* 

- 


_ j 

- 

186.5 


- 

- 

- 

! 356.8 


345.3* 


i 

“ 1 

- 

- 

227.9 


- 

- 

' 400 

1 

— 

— 

422.0* 

478.9* 

i 

518.7* 

~ 

— 

288.0 

328.3 

356.7 


*Floal Rise 
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TABLE 8-6 


PLI ME RISE PARAMETERS FOR M(H)EBATE WIND GONDITICSS 



i 

♦ 

Centerline Height (a) 


Pluae Radios (a) 

DoHnwind: 

Distance 

Huaber of Sega«its 


Huaber of Segaents 

(*) 


2 

4 

1 

6 ; 

-1 

8 

D 

2 1 


■1 

8 

100 

72.1 

93.0 

103.4 

1 

114.9 

122.9 

51.3 

61.4 ; 

77.7 

88.1 

95.5 

125 

82.8 

105.9 

119.4 

133.1 

142.5 

58.3 

69.9 

88.3 

100.1 


160 

|%.3 

122.2 

139.7 


167.3 

67.2 

80.7 


115.2 

124.7 


110.2 

139.0 

160.6 

179.7 

192.9 

76.4 

91.8 

115.4 

130.8 

141.7 


Mi 

158.0 

184.2 

206.5 

221.9 

86.8 

104.3 

131.1 

148.5 



145.7 

181.9 

213.9 


258.2 

99.8 

120.0 

150.7 


184.7 

1 400 

165.8 

206.2 

244.2 

274.4 

295.4 

113.1 

136.1 

170.7 

193.3 



188,1 

233.2 

278.0 

312.6 

336.8 

127.8 

153.9 

192.9 

218.6 

236.5 


213.4 

263.9 

316.4 

356.1 

384.1 

144.5 

174.2 

218.3 

247.4 

267.9 


241.4 

298.1 

359.4 

405.1 

437.1 

162.9 

196.8 

246.6 

279.6 

302.7 


268.0 

331.0 

401.0 

452.7 

488.8 

180.6 

218.5 

274.1 

311.0 

336.9 


293.1 

362.4 

441.2 

499.0 

539.5 

197.1 

239.2 

300.6 

341.6 

370.4 


313.3 

389.1 

477.1 

541.4 

586.4 

210.5 

256.8 

324.4 

369.5 


1828.7 

317.2* 

- 

- 

- 

- 



- 

- 

- 

1889.6 


396.2* 

- 

- 

- 

- 

261.5 

- 

- 

- 

1987.7 

- 

- 

- 

557.9* 

- 

- 

- 

- 

380.5 

- 

■ 

— 

— 

489.4* 

1 

— 

606.2* 


— 

332.4 

1 

414.5 


♦Final Rise 
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(Equation (9-8) continued) 



where 


Q 


source emission rate 


a 

y 


standard deviation of the crosswind concentration 
distribution 




ry 


/ X + X - X (1-a) 

!: ^ — 

\ ry 


( 8 - 9 ) 


oMt} 

A 

o!{t} 

A 


T 


standard deviation of the wind azimuth angle in radians 
for the averaging time T 


(I 




( 8 - 10 ) 


"a • 

T • T > T 
E ’ A E 


( 8 - 11 ) 


A 0 


X 


A 


X 


standard deviation of the wind azimuth angle in radians 
for the meteorological sampling time period x^ 

concentration averaging time 

downwind distance from the source 
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crossvii^ vircual distance 



crosswind dispersion coefficient 


( 8 - 12 ) 


distance over which rectilinear crosswind cloud expansion 
occiirs downvind froa a virtual point source 


standard deviation of the crosswind concentration distri- 
bution at a distance x_ fro^ the source 


standard deviation of the vertical concentration 
distribution 


a- (, + x^) 


( 8 - 13 ) 


standard deviation of the wind elevation angle in 
radians 


vertical virtual distancv 


zR 

o: ■ ’Si 


( 8 - 14 ) 


standard deviation of the vertical concentration 
distribution at 


laean wind speed 

1+p l+p 

^ -K) 

h ■ "i) isY 


; u > u. 


; u < u 


( 8 - 15 ) 


93 



aean wind speed ac the reference height 


wind profile power-law exponent 


effective lower bound of the cloud 


H + 2.15 o ; z, > 2 
z I 


; * 1^2 


effective upper bound of the clcaid 


H + 2.15 a ; z, < H 
z 2 o 


= '2 i 


crosswind distance from the cloud centerline 


depth of the surface mixing layer 


effective height of the source above ground level 


Ah from Equation (8-1) or (8-5) 


height above ground level 


time after release 


standard deviation of the alongwind c •>ncencratlon 
distribuClcn 
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1/2 


a 

X 



2 

o 

xo 


L{x} 


0.6 Au 


x ; lu i. 0 


u < 0 


( 8 - 19 ) 


(8-20) 


Au = vertical wind-speed shear ir the layer containing the 
cloud 



( 8 - 21 ) 


a = 
xo 


standard deviation of the alongwind cloud distribution 
at the source 



( 8 - 22 ) 


( 8 - 23 ) 
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Calculation Procedures 


The concentrations for downwind distances between the source and 
the distances of maxlmua plume rise shown in Table 8-S and 8-6 for the 
various segment configurations were calculated using Equation (8-8) at each 
distance x using the meteorological input parameters in Table 8-3, the 
source parameters in Table 8-4, and the following additional input parameters: 



The values of o , a , a and at the distance of maximum 
yr zr’ xo K 

plume rise were used to calculate concentrations using Equation (8-8) for all 
distances be>ond the point of maximum cloud rise. The concentration calcula- 
tions were made for nearly-instaneous peak concentrations ™ 2.5s) and 

A 

for lO minute (,T = 600s) and 30 minute (r = 1800s) averaging times. 

A A 
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Results of the Calculations 


The results of the calculations are shown in Figures 8-1 through 8-10. 
Figure 8-1 shows Instantaneous and 10- and 30-niinute time average HCl concen- 
trations downwind from a l-segment bum during very light wind conditions. 

At distances close to the source. Figure 8-1 shows that the Instantaneous 
ana 10-minute time average concentrations are identical, a result that we 
expect because the segment burns for about 15.5 minutes ac a constant oolssion 
rate. The 30-minuto time average concentration is less because concentra- 
tions are nearly zero during part of the 30 minutes after the terminated 
plume has passed the receptor. At these distances close to the source, 
concentration decreases with increasing distance because the plume is rising 
faster than it is expanding and the plume begins to lose contact with the 
ground. As turbulence continues to cause plume expansion after maximum 
plume rise has been achieved, concentrations again begin tu increase to a 
maximum value. At longer downwind distances when the terminated plume is 
expanding alongwind due to wind speed shear in the vertical, the 10-mlnute 
and 30-minute time averaged concentrations become more nearly identical 
(the plume is "long" compared to the averaging time). The remaining figures 
for burns of the various segment configurations show similar patterns. 

The highest concentrations beyond 1 kilometer from the sources occur under 
"worst-case" (very light wind conditions) and are greatest for the 8-segment 
(2SRM's) storage configuration reaching a maximum instantaneous concentration 
of 66 ppm HCl at 6.5 kilometers downwind from the storage area. 
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FIGURE 8-4. Instantaneous HCl concentrations and 10- and 30-ninuts tins average HCl concentrations 
downwind from a 2-segment burn during moderate wind conditions. 
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FIGURE 6-5. Inscantaneous HCl concencratlons and 10- and 30-vinute Clna avtrage ilCl coocancratlona 
dotmwlnd from a 4-aagmant (ISRN) burn during very light wind conditions. 
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FIGUKE 8-6. Instantaneous HCl concentrations and 10- and 30-nlnute tiaa average HCl concentretlone 
doimwind from a 4-segnent (ISRH) burn during moderate wind conditions. 




FIGURE 8-7. Instantaneous HCI concentrations and 10- and 30-nlnut« tlisa average HCl concentrations 
downwind from a 6-segment bum during very light wind conditions. 




$ 


FIGURE 8-8. Instantaneous HCl concentrations and 10- and 30-mlnute tine average concentrations 
downwind from a 6-segment bum during moderate wind conditions. 












(ludd) N0liV«lN30N0D IDH 
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FIGURE 8-9. Instantaneous HOI concentrations and 10- and 30-nlnute time average HCl concent vat ions 
dovmwlnd from a 8-scgricnt (ZSRll'u) bunt during very light wind conditions. 
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FIGURE 8-10. Instantaneous HCl concentrations and 10- and30-minute time average HCl concentrations 
downwind from a 8-sef.t ent (2SRM's) burn during moderate wind conditions. 



SECnON 9 


HAZARD CALCULATI(»S m POSTOLAIED SPILLS OF MMI-10 
FDEL AT FOQ. FARM fl, KSC 


Instantaneous and tlae-av^rage NO 2 concentratioR profiles have 
beat calculated downwind froa postulated spills of MOR-IO fuel at Fuel 
Fan #1 at KSC. The calculations tiere perforated for three Beteorological 
cases specified by DF-EM5/ISC. In Case 1, the wind at a height of 2 a is 
assist to be froa 199° at a speed of 3 a s According to KSC, this 
wind direction takes tt» evaporating cloud towards the CCAFS Industrial 
Area at a distance of 3353 a froa Fuel Fan #1. In Case 2 , a wind of 4 
a s ^ fm 130° is asstaed to carry the cloiki towards the KSC Industrial 
Area located at a distance of 8534 a froa Fi»l Fan #1. In Case 3, a 
wind froa 019° at a speed of 5 a s ^ is assined to carry the cloud towards 
Cape Canaveral Harbor at a distance of 6706 aeters froa Fuel Fan #1. 

In each case, DF-EIC/KSC suggested chat we assume that the MON-10 

2 2 

fuel covers an area of 1600 ft (148.65 a ) within Che dikes constructed 
at Fuel Farm #1. Recognizing that MON- 10 fuel is comprised of 90Z ^20^ and 
lOZ NO, it was also suggest^ chat a spill of 27,675 pounds of ^20^ and 
3,075 pounds of tK) be considered for use in Che calculations. After dis- 
cussion with DF-EMS/KSC, a joint decision was made to consider only the 
conversion of the 27,675 pounds of N20^ Co NO 2 gas in the calculations. 

Also, DF-EMS/KSC suggested an evaporation rate for MON-10 fuel of 5 x 10 ^ 

gal min ^ ft Since MON-10 fuel has a density of 12.3 lbs gal ^ at 68° F, 

-2 -1 

this is equivalent to an evaporation rate of 5.0045 g m s . 


Meteorological Model Innut Parameters 

As noted above, the wind conditions for the calculations were 
specified by DF-EMS/KSC. The other meteorological input parameters we 
selected as being representative of meteorological conditions at KSC under 
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tlw specified wind conditions are shown in Table 9>1. The yaliMs of the 
wind power-l«f exponent p and the turbulence paraaeters and in 
Tabic 9-1 are conservative values based on an analysis of aeteorological 
data fnm the NASA-150 seter tower at KSC (see Kecord, et 1969). That 
ls» the use of these values In the aodel calculations should lead to pre- 
dicted concentrations that are aavliw values for the specified wind con- 
ditions. The values of tlw surface nixing depth in Table 9-1 are 
considered to be average values for these wind conditions at KSC. There 
nay be tlaes» particularly late at night and early in the nomii^ with 
winds of 3 a s ^ and less, when the mixing depth is lower and would act 
to increase co.-'centrations at longer travel distances. The values of 
the dispersion oiefficients a and 6 given in Table 9-1 are based on 
our experience in modeling concentrations downwind from quasi-continuous 
sources. 


Source Model Input Parameters 


Tl« source model input parameters used in the calculations are 

given in Table 9-2. The values of the initial lateral (o _} and alongwind 

yK 2 

((7^^) source diaeaslons were obtained by assuming the spill area of 148.65 m 
to be circular and dividing the diameter of 13.7574 a by a factor of 4.3 to 
obtain the standard deviation of the concoitration distribution at the 
source. Thus, 



13.7574 

4.3 


3.2 a 


The vertical source dimension was arbitrarily set to 0.1 a. The source 

emission time T_ is 
c* 


T 


E 


27,675 lbs 


X 


453. 6r b^s 

a X 

lb 5.0045 g 


1 

148.65 m^ 


- 16,874 s 
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TABLE 9-1 

METEOROLOGICAL MODEL INPUT PARAMETERS 


Model Parameters 

Case 1 



Case 2 

Case 3 

u {2m} (m s 

3 

4 

5 

P 

0.175 

0.15 

0.12 

H Ob) 

m 

600 

800 

800 

0. (t *10 min} (d^) 

A O 

8 

8 

8 

Og (deg) 

2.7 

2.7 

2.7 

a 

0.9 

0.9 

0.9 

6 

1.0 

1.0 

1.0 
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TABLE 9>2 

SOURCE MCfflEL INPUT PARAMETERS 


Model Paraaeters 

Value 

V ■ “k, 

3.2 

®zR 

0.1 


0 

V - <■> 

50 

Tj Cs) 

16,874 

Q (ppm m^) 

6.6708 X 10^ 

H (m) 

0 
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The amount of NO^ that can be obtained from 27,675 pounds of N20^ is 


Q {NO2} 


27,675 lbs N20^ 


X 

Ab.OlgNO^ 


m^ 298.16°K 453. 6 r 

— .1 II J 

lo\ 273.16°K lb 


6.6708 X 


m3 3 
10 m 


Since we wish the concentration to be given in parts of NO. per million parts 
of air, the above value of Q must be multiplied by 10 , or 

Q {Q} * 6.6708 X 10^ X 10^ * 6.6708 x 10^ ppm m^ 

which is the value for Q shown in Table 9-2. Finally, the source release 
height H is set to zero for a ground-level release. The calculations were 
perfoimed for nearly- instantaneous (T^ = 2.5s), 10-minute (T^ = 600s) and 
30-mlnute (T^ = 1800s) time-average concentrations. 

Results of the Calculations 


The concentration calculations were made using the quasl-continuous 
source model described above in Section 8. Figures 9-1 through 9-3 show 
calculated peak NO 2 concentrations versus downwind distance from the hON-10 
fuel spill for Cases 1 through 3. In each figure, the three concentration 
curves represent averaging times of 2.5 seconds (Instantaneous), 10- and 
30-mlnutes. The NO 2 concentrations calculated at the distances of the 
CCAFS Industrial a, ia (Case 1), KSC Industrial Area (Case 2), Cape Canaveral 
Harbor (Case 3) are presented in Table 9-3. The calculated lO-mlnute and 
30-minute NO 2 concentrations at the KSC Industrial Area (Case 2) are below 
the lowest concentration shown in Figure 9-2. 

Figure 9-4, 9-5 and 9-6 show the time profiles of instantaneous 
concentration for Cases 1 through 3 at the CCAFS Industrial Area, the KSC 
Industrial Area and Cape Canaveral Harbor. Figure 4 shows, for example. 
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NO* CONCENTRATION (ppm) 



FIGURE 9-1. Instantaneous, 10- and SO-aicute tiae average NO 2 ccncentrations 
for a 2-neter wind speed of 3 meters per second (Case 1). 
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UOz CONCENTRATION (ppm) 
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NOj CONCENTRATION (ppm) 



DOWNWIND DISTANCE (km) 


FIGURZ 9-3. Instantaneous, 10- and 30-nlriuto tloc average NO 2 concentrations 
for a 2-«eter wind speed of 5 meters per second (Case 3). 
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TABLE 9>3 


CALCULATED NO2 COHCENTRAIIOMS AT THE CCAFS INDCSTRLAL AREA, 
KSC IMDUSTRIAL AREA AMD AX CAPE CAMAVERAL HARBMl 


i 

NO2 Concentration (ppa) 

1 

Instantaneous 

10-ain 

30-ain 

CCAFS Industrial Area 
(Case 1, X « 3353a) 

l.l 

0.37 

0.30 

KSC Industrial Area 
(Case 2f X * 8534a) 

0.14 

O.Oj 

0.04 

1 Cape Canaveral Harbor 
(Case 3. X * 6706 a) 

i 

0.21 

0.07 

0.06 
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FlCUiU; 9-4, Time profile of instantaneous NO 2 concentrations at the CCAFS 

Industrial Area for a 2-meter wind speed of 3 icaters per second 
(Case 1). 
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FIGURE 9-5. Tine profile of Instantaneous NO 2 concentrations at t^a SaC 

Industrial Area for a 2-swter vind speed of 4 meters per second 
(Case 2). 



CONCENTRATtON {ppm) 












that the concentration exceeds 1 ppa at the CCAFS Industrial Area about 
10 minutes after the spill occurs and remains greater, on the average , than 
1 ppa until 289 minutes after the spill occurs. The time of arrival and 
the passage time of various Instantaneous concentration levels at the 
CCAFS and KSC Industrial Areas and for Cape Canaveral Harbor can be 
estimated frcmi the curves In the three figures. 
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